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Abstract 


Experiments have been conducted to investigate the 
effect of axial jet-induced mixing on the pressure reduc- 
tion of a thermally stratified liquid hydrogen tank. The 
tank was nearly cylindrical, having a volume of about 
0.144 m 3 with 0.559 m in diameter and 0.711 m long. A 
mixer/pump unit, which had a jet nozzle outlet of 
0.0221 m in diameter was located 0.178 m from the tank 
bottom and was installed inside the tank to generate the 
axial jet mixing and tank fluid circulation. The liquid 
fill and jet flow rate ranged from 42 to 85 percent (by 
volume) and 0.409 to 2.43 m 3 /hr, respectively. Mixing 
tests began with the tank pressure ranging from 187.5 to 
238.5 kPa at which the thermal stratification results in 
4.9 to 6.2 K liquid subcooling. The mixing time and 
transient vapor condensation rate at the liquid- vapor 
interface are determined. Two mixing time correlations, 
based on the thermal equilibrium and pressure equilib- 
rium, are developed. Both mixing time correlations are 
expressed as functions of system and buoyancy param- 
eters and compared well with other experimental data. 
The steady state condensation rate correlation of Sonin 
et al. based on steam-water data is modified and ex- 
pressed as a function of jet subcooling. The limited liq- 
uid hydrogen data of the present study shows that the 
modified steady state condensation rate correlation may 
be used to predict the transient condensation rate in a 
mxing process if the instantaneous values of the jet sub- 
cooling and the turbulence intensity at the interface are 
employed. 


Nomenclature 


A 8 liquid-vapor interface area 

C p specific heat of liquid 

D tank diameter 

d jet nozzle diameter 

g gravitational acceleration 

H liquid height from the tank bottom 

Hj jet submergence depth or the distance be- 

tween nozzle outlet and liquid surface 


height measured from the tank bottom 
latent heat of condensation 
liquid Jakob number, C p (T g - T b )/h fg 
tank height 

average mass condensation rate at the 
liquid-vapor interface 

tank pressusre 

Prandtl number of bulk liquid 
jet volume flow rate 

Richardson number evaluated at the be- 
ginning of each test, g/3(T B — Tj)j Hj/uj 2 

jet Reynolds number, Uj d/i/ 

average condensation Stanton number de- 
fined in Eq. (8) 

bulk liquid temperature, jet temperature, 
interace temperature, and dome 
temperature 

initial jet subcooling of each test 

mixing time, dimensionless mixing time 

mixing time based on the thermal and 
pressure equilibrium, respectively 

jet velocity 

liquid volume 

v g average rms turbulent velocity (or turbu- 

lence intensity) at the liquid-vapor 
interface 

Greek Symbols 

P thermal expansion of liquid 

P V P 2 constants defined in Eq. (10) 

v kinetic viscosity of liquid 

p liquid density 
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Introduction 


Okit& and Oyama 3 


Fluid mixing such as axial jet-induced mixing is a 
preferred method to reduce the tank pressure of a two- 
phase system in a closed container in which the fluid is 
thermally stratified. 1 A typical example is the thermal 
destratification and the pressure reduction of a cryogenic 
storage tank by axial jet-induced mixing. As the tank is 
subjected to external heating, the fluid is thermally stra- 
tified and the tank pressure increases with a rate usually 
higher than its homogeneous pressure rise rate. By using 
a mixer, the tank fluid is circulated and the subcooled 
jet can be directed towards the liquid-vapor interface. 
The mixing of the tank fluid induced by the jet can de- 
stroy the temperature stratification and induce the 
vapor condensation at the interface, resulting in the 
reduction of tank pressure. Once the liquid temperature 
stratification is eliminated, the bulk liquid temperature 
is nearly saturated. Consequently, the continuation of 
mixing will not further reduce the tank pressure. In- 
stead, the mixer power input, which may be significant 
for a cryogenic storage system, will eventually become 
heat and increase the boiloff losses. Therefore, the mix- 
ing process, which involves the mixing time and the va- 
por condensation rate at the interface, needs to be well 
understood such that the pressure control system can be 
properly designed and operated. 

Fluid mixing process has been extensively investi- 
gated and numerous correlations to predict the mixing 
time have been proposed. However, the concept to for- 
mulate the mixing time data differ from one to another. 
Most of the tests were conducted with cylindrical tanks 
using water as the test fluid with a uniform temperature 
and injected dye of equal density to eliminate the buoy- 
ancy force. Some of the previous studies on mixing time 
are as follows: 
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For comparison purpose, Okita and Oyama 3 recalculated 
the data of Fox and Gex 2 and obtained the follow- 
ing correlation which is independent of jet Reynolds 
number: 
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for Rej > 7000 ( 5 ) 


From the above correlations, it is found that for 
H/D close to one the dimensionless mixing time, t m 
== t m u. d/D 2 , ranged from 2.6 to 9. Since all the above 
tests used the dye or acid-base method, the mixing time 
results were based on the concentration equilibrium 
which may be comparable to the thermal equilibrium. 

Poth and Van Hook 6 conducted a mixing test by 
using liquid hydrogen stored in the Nuclear Propulsion 
Module (NPM) Mars Braking Stage. The mixing time 
for the tank pressure to reach its equilibrium value is 
about 
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Similar experiments were conducted by Meserole et al. 
using Freon 11 at various jet diameters and liquid 
heights. They concluded that the mixing time deter- 
mined from the pressure decay is essentially equal to 
that based on the equalization of temperatures through- 
out the bulk liquid. Also, the mixing time only varies 
slightly with fill levels. It seems that their data may be 
correlated as 


i„=ii85^_i L 

ReV 6 g l/. 


for Rej > 2000 


( 2 ) 


t = 3.9 £_ for Re; > 20 000 (7) 

® Ujd J 


2 



The reduction of tank pressure is mainly due to the 
vapor condensation at the interface. The mixing time 
definition based on the pressure equilibrium should be 
affected by the pressure decay rate which is controlled 
by the vapor condensation rate at the interface and the 
thermodynamic state of the vapor. 

Several experiments have been conducted for the 
measurement of the steady state vapor condensation rate 
at the liquid-vapor interface in a jet-induced mixing 
tank. Thomas 7 measured the condensation rate of steam 
on water surfaces mixed by a submerged turbulent jet 
and attempted to correlate his data based on the theo- 
retical model proposed by Theofanous et al®. Theofanous 
model suggests that the condensation heat transfer 
coefficient be proportional to the Reynolds number to a 
power of 0.75. Dominick 9 investigated the effects of jet 
injection angle and jet flow rate on the condensation 
rate in a Freon 113 tank. The interface heat transfer was 
observed to increase as jet injection angle became more 
normal to the interface. Dominick 9 correlated his data 
by using the correlation form which was developed for 
a turbulent jet impinging on a flat surface. In this 
correlation, the average heat transfer coefficient increases 
with increasing jet Reynolds number to a power of 0.73. 
Sonin et al. 10 obtained an empirical correlation for the 
rate at which a pure vapor condenses on the free surface 
of a turbulent liquid in a steam-water system, for the 
conditions where buoyancy effects are insignificant. The 
vapor was saturated and the vapor condensation was 
controlled by the interface turbulence at the liquid side. 
Their vapor condensation rate correlation is expressed as 

St = m c h fg 

V C p v s( T . - T b) (8) 

= 0.0198 Pr -0 33 (1 - 0.5 Ja) 

where St c is the average condensation Stanton number 
and m c is the average vapor mass condensation rate at 
the interface. The liquid subcooling (T g - T b ) and 
average turbulence intensity at the interface (v g ) are 
mainly responsible for the vapor condensation rate. The 
average turbulent rms velocity at the interface, v , is 
determined by 
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where = 0.34 and = 0.24. Lin and Hasan*** 
numerically solved the conservation equations with k-e 
turbulence model and obtained that 


v = ^ 2.39 - 0.645 for 0.5 < 5 < 3 ( U ) 

D ( Dj D 

Equation (11) yields = 0.33 and fl 2 = 0*23 which are 
in good agreement with the results of Ref. 12. It is noted 
that Eqs. (9) to (11) were obtained for a cylindrical tank 
with the jet nozzle outlet located near the tank bottom. 
The liquid height from the tank bottom (H) is nearly 
the same as jet submergence depth (H.). For most of the 
practical applications, they are not the same and Hj is 
suggested to be used instead of H. 

This paper presents the experimental results of the 
pressure decay of a liquid hydrogen tank mixed by an 
axial turbulent jet directed normally towards the liquid- 
vapor interface. The controlling variables are the jet 
volume flow rate and liquid fill level. Since the tank 
pressure is decreasing during the fluid mixing, the con- 
densation rate is changing with time because of the 
change in jet subcooling. The steady state condensation 
rate correlation developed by Sonin et al. 10 is modified 
to be expressed as a function of jet subcooling (T g - Tj). 
The transient condensation rate during the mixing 
process is then predicted by this modified steady state 
condsensation rate correlation using the instantaneous 
jet subcooling. The correlations of the mixing time, 
based on both the thermal equilibrium and the pressure 
equilibrium, are developed and are compared with other 
experimental data. 


The Experiment 

The mixing tests are performed in a nearly cylin- 
drical dewar. The dewar is a vacuum jacketed stainless 
steel tank containing multilayer insulation (MLI) within 
the vacuum annulus. The lid is composed of a flat flange 
which supports a short cylindrical section with an 
inverted dome bottom. The space between the flange 
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and the cylindrical section is evacuated and insulated 
with MLI to minimize heat leak through the dome from 
the environment. With the lid in place, the interior walls 
of assembled tank form a cylindrical storage volume 
with dome ends. The assembled tank then has an inter- 
nal volume of about 0.144 m 3 with an internal height of 
0.711 m and an inside diameter of 0.559 m. Because of 
the lid mounted penetrations and the coupling of the lid 
walls to ambient temperatures at the tank flange, the 
heat leak is very nonuniform. The average wall heat flux 
for the tank is experimentally determined and ranges 
from 3.15 to 31.5 W/m 2 for fill levels ranging from 30 to 
85 percent. The Shuttle Centaur mixer is installed near 
the central part of the tank bottom with the jet nozzle 
outlet located about 0.178 m from the bottiom. The jet 
nozzle has an inner diameter of 0.0221 m. With the 
mixer on, the tank fluid is withdrawn and ejected 
through the jet nozzle normally towards the liquid-vapor 
interface, resulting in axial jet-induced mixing in the 
tank. Figure 1(a) shows a rough sketch of the tank. It is 
noted that the dewar used in the present mixing test is 
indeed the receiver dewar of the no-vent fill tests in 
Ref. 14. The details of the structure, piping systems and 
the data acquisition system can be found in Ref. 14. 

Liquid height is measured by a capacitance level 
probe. An instrument tree containing silicon diodes and 
thermistors at varying heights is in direct contact with 
the tank fluid. The accuracy of the silicon diodes is 
about 0.6 K at LH 2 temperatures. Six silicon diodes are 
located at 0.5-, 20-, 40-, 50-, 60-, and 80-percent fill 
height to measure the fluid temperatures. Four thermis- 
tors are utilized to indicate the presence of liquid or 
vapor. The capacitance level probe is calibrated against 
the thermistors (point sensors) and is found to agree 
within 0.0254 m for liquid levels greater than 10 percent. 
There are two additional silicon diodes used to measure 
the tank dome temperatures. Figure 1(b) illustrates the 
temperature sensor and thermistor locations for the test 
tank. Pressure transducers are used to provide continu- 
ous pressure measurement throughout the system with 
an estimated accuracy of ±1.7 kPa. A silicon diode and 
a flow meter are also installed in the transfer line of the 
mixer unit to measure the jet temperature and jet 
volume flow rate, respectively. 

Before each test, the tank is vented to 40 to 60 kPa 
at a specified liquid fill level. In order to establish the 
initial conditions for each test, the tank vent valve is 
then closed. Due to the ambient heat flux through the 
tank wall, the tank pressure increases and the liquid 
becomes thermally stratified. When the tank pressure 
reaches a preset level, the mixer is turned on at a speci- 
fied jet volume flow rate and the test is considered to 
begin. The data for the liquid fill level, fluid tempera- 
tures, tank pressure, and jet flow rate are measured 
starting just prior to activating the mixer. During each 


test, the sampling rates of the measurement are varied 
from once every 2 sec to once every 10 sec. The mixer is 
turned off when the pressure measurement does not 
show any pressure decrease for a specified period of time. 
The vent valve is then open and the test is terminated. 

Results and Discussion 

Four series of experiments consisting of 12 tests have 
been conducted. Table 1 shows the experimental condi- 
tions for all the tests. The liquid fill ranged from 42 to 
85 percent (by volume) and the initial liquid subcooling 
due to thermal stratification ranged from 4.9 to 6.2 K. 
The jet flow rates varied from 0.409 to 2.43 m 3 /hr. 
During each test period, the liquid fill level and the 
dome temperature hardly changed. Therefore, only their 
initial values are shown in the table. 

Figures 2(a) to (d) show the tank pressure decay as 
a function of time for various jet volume flow rates. It 
is seen that the pressure decay rate is higher for the 
higher jet flow rate as the result of the enhanced vapor 
condensation at the liquid-vapor interface. Also note in 
Figs. 2(a) to (d), the lowest jet flow rates substantially 
increase the mixing time. The typical liquid temperature 
distributions as functions of time are shown in Figs. 3(a) 
and (b). The liquid temperature is completely destrati- 
fied and the tank pressure reaches its equilibrium value. 
This equilibrium pressure is nearly the same as the 
saturation pressure corresponding to the uniform bulk 
liquid temperature of the tank. For the tests of the same 
series, the tank almost reaches the same equilibrium 
pressure. This is because their initial experimental con- 
ditions such as liquid fill, tank pressure, and fluid tem- 
perature distributions are nearly the same and thus 
should have the same equilibrium state when the tank 
is completely mixed. 

The mixing characteristics for all the tests are given 
in Table 2. The jet Reynolds number ranges from 27 900 
to 166 300. A Richardson number, Ri, is also included 
in Table 2 to characterize the effect of buoyancy. The 
definition of the Richardson number is based on the ini- 
tial jet subcooling (T # - Tj)., jet submergence depth 
(H-), and jet velocity (u.). The value of Richardson 
number ranges from 0.07 to 2.62. Two mixing times are 
shown in the table: (t m ) th and (t m ) p . The mixing time, 
(t m ) th , is defined as the time required from the start of 
a mixing operation to the time when the liquid tempera- 
ture is destratified. The mixing time, (t m ) p , is defined as 
the time required from the start of a mixing operation 
to the time when the equilibrium tank pressure is 
achieved. Table 2 shows that, for small Ri, the dimen- 
sionless mixing time (t m ) th is independent of the^jet 
Reynolds number. The dimensionless mixing time (t m ) th 
as a function of Richardson number is shown in Fig. 4. 
As shown in Fig. 4 the effect of buoyancy on the mixing 
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time is negligible for Ri < 0.5. The small value of Ri 
can be obtained by increasing the jet velocity (uj) or 
decreasing the gravity level (g). Therefore, the dimen- 
sionless mixing time for the high jet flow rates in a 
normal-gravity condition may be applied to the low jet 
flow rates in a zero-gravity environment, assuming the 
effect of liquid-vapor interface configuration is neglected. 
Using the form of the correlation developed by Okita 
and Oyama 3 , the mixing time data shown in Table 2 
can be correlated as 



(‘■)th«J d 

D 2 



Ri < 0.5 (12a) 


(3.2 + 2.0 Ri) 


'«r 




Ri > 0.5 


(12b) 


may be different as shown in Table 2. If the tank is 
completely mixed and the bulk liquid is still subcooled, 
the condensation process and the pressure decrease 
continue. In this case, (t m ) th is less than (t m ) p . The low 
liquid fill data shows this phenomena. If the condensa- 
tion process stops before the tank is completely mixed, 
then (t m ) th is greater than (t m ) p . This may happen for 
very high jet volume flow rates. It is obvious that the 
value of (t m ) p can not be correlated simply by the sys- 
tem parameters as those expressed in Eq. (12). The 
value of (t m ) p depends on the pressure decay rate which 
is controlled by the condensation rate at the liquid- 
vapor interface and the thermodynamic state of the 
vapor region (particularly the vapor mass). This can be 
examined by simply looking at the equation of the 
pressure change rate for an isothermal vapor region. 
However, for comparison purpose, the value of (t m ) p as 
a function of Ri is plotted in Fig. 5. Similar to that for 
(t m ) th , buoyancy increases the value of (t m ) p . With neg- 
ligible buoyancy, the value of (t m ) p increases with 
decreasing fill level. This is because the larger vapor 
mass for the low fill level delays the pressure decay rate. 
This phenomena was also observed by Meserole et al. 
but was not discussed. Our data may be correlated by 
the following equation: 


Equation (12a) shows that, for Ri < 0.5, the mixing 
time (t m ) th decreases linearly with increasing jet velocity 
for a given tank system and decreases with decreasing 
fill level. This finding is consistent with that of Okita 
and Oyama. 3 However, the constant 4.3 in Eq. (12a) for 
our liquid hydrogen system is smaller than the con- 
stant 5.5 in Eq. (3) obtained by Okita and Oyama 3 for 
a water system. For higher Richardson number 
(Ri > 0.5), the mixing time is higher (as shown in 
Fig. 4) due to the opposing effect of the buoyancy. 
Equation (12b) gives the quantitative effect of buoyancy 
on the mixing time. It is noted that in terms of jet 
volume flow rate and total liquid volume, Eqs. (12a) 
and (12b) can be rearranged to be of the form: 


Ri < 0.5 (13a) 


Ri > 0.5 (13b) 


Since the tank pressure is controlled by what hap- 
pens in the vapor region, the values of (t m ) th and (t m ) p 
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For Ri less than 0.5, the constant 3.3 in Eq. (14a) is 
slightly less than the constant 3.9 obtained from 
Meserole et al.'s 7 data for a Freon 11 system. It is 
noted that for a space system in which the gravity effect 
is neglected, Eqs. (12a) and (14a) may be used to 
predict the mixing times (t m ) th and (t m ) p , respectively. 

Vapor condensation at the interface is the key 
mechanism to reduce the tank pressure. For a well 
insulated tank, the interface condensation processes are 
generally determined by the system parameters (such as 
tank geometry and size, jet nozzle size, and liquid 
height), fluid properties, liquid subcooling, jet flow rate, 
and buoyancy. Figures 6(a) and (b) show the jet sub- 
cooling as a function of time for high jet flow rates. 
With the corresponding change of tank pressure shown 
in Figs. 2(b) and (d), two conclusions can be drawn. 
First, for a given jet flow rate, pressure decay rate 
decreases with time due to the decrease in jet subcool- 
ing. Secondly, for a given jet subcooling, higher jet flow 
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yields a greater pressure decay rate. This is because the 
vapor condensation rate increases with increasing jet 
subcooling and with the enhanced turbulence intensity 
near the interface (due to the increased jet flow rate). 

The vapor condensatin rate correlation developed by 
Sonin et al. 10 is expressed as 


Using Sonin et al.'s 10 correlation, St c — 0.0198 Pr °* 33 
we have 
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m c h fg 
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= 0.0198 Pr _0,33 (l - 0.5 Ja) 


Since this correlation is for steady state, there is a direct 
relationship between the interface temperature (T g ), 
bulk liquid temperature (Tj, and jet temperature (Tj). 
The subcooled jet is the real driver to induce the vapor 
coondensation. The subcooled jet (Tj) approaches and 
sweeps over the interface, and induces the vapor conden- 
sation at the interface. The jet flow absorbs the energy 
through this condensation process and leaves the inter- 
face at a temperature T b . Therefore, the bulk liquid 
temperature (T b ) is greater than the jet temperature 
(Tj) and is less than the interface temperature (T g ). In 
the mixing test, the condensation process is transient be- 
cause the interface temperature keeps decreasing. It may 
be difficult to measure the temperature T b because the 
condensation-induced thermal layer is growing and the 
temperature T b is changing with time. Therefore, it will 
be helpful if the liquid subcooling (T a - T b ) in Eq. (15) 
can be expressed as a function of the jet subcooling (T g 
— Tj) which can be easily obtained in the experiment. 

Based on the mass and energy balance, Lin and 
Hasan 16 have theoretically obtained the following rela- 
tionship between (T 8 — Tj) and (T b — Tj) for an insu- 
lated system: 
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Equation (18) indicates that the diameter ratio (D/d) is 
the key parameter controlling the difference between the 
liquid subcooling(T 4 - T b ) and jet subcooling (T, 
— T-). For the system used in the present experiment, 
the value of (T a — T b ) is nearly half of (T a — Ts), if the 
value of v a is calculated by Eq. (11)* Using Eq. (18), 
Sonin et al.'s 10 condensation rate correlation becomes 


st c = 


m c hfg 


A 8 pC p v a (T a - Tj) 

(0.0198 p r ~ 0,33 ) 
= 0.0198 Pr~°* 33 (l - 0.5 Ja) 
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D 


J. 


(19) 


where the fluid properties are calculated based on the 
temperature T b which is obtained from Eq. (18). It is 
noted that Eq. (18) holds for a steady state condition. 
Assuming Sonin et al.'s 10 condensation rate correlation 
is applicable to a transient state, the instantaneous val- 
ues of (T a - T b ), (T a — Tj), and v a should be used. Un- 
fortunately, the information to predict v 8 as a function 
of time is not available. Theoretically, the steady state 
value of v a can be used only when the tank has been 
completely mixed and the required mixing time is that 
predicted by Eq. (12a). This steady state value of v g can 
be determined by Eqs. (9) and (11) for high H/D and 
low H/D, respectively. The average turbulence intensity 
(v g ) is expected to increase with time in a mixing 
process. If one uses the steady state value of v g , the 
condensation rate is expected to be overpredicted at 
least in the initial time period of mixing process. This 
effect of transient v g may be decreased with increasing 
jet flow rate because the mixing time is shorter. 

Based on the limited temperature measurements in 
the vapor region, curve fitting is used to estimate the 
vapor temperature distribution as a function of time. 
The typical fluid temperature distribution as a function 
of time is shown in Figs. 7(a) and (b). The vapor is su- 
perheated with the temperature difference (T dome - T g ) 
equal to about 65 K. The solid line is the resulting 
curve-fitted temperature distributions at the time when 
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the mixing process starts. The vapor mass is obtained 
based on the curve-fitted one-dimentional temperature 
distribution. Most of the vapor mass is confined in the 
region near the interface where the vapor temperatures 
are much lower. Figures 8(a) to (e) show the estimated 
vapor mass as a function of time for the high jet flow 
rates (Ri < 0.5). The vapor mass predicted by using 
Eq. (19) with the steady state value of v # is also in- 
cluded in the figures for comparison. As expected the 
condensation rate is generally overpredicted. This is 
mainly due to the use of the steady state value of v g and 
slightly due to the vapor superheat. The comparison ap- 
pears to improve as the jet flow rate increases. It is 
noted that the condensation rate is underpredicted at 
the jet flow rate of 2.43 m 3 * /hr as shown in Fig. 8(c), 
which is the highest flow rate used in the experiment. 
This may be because at this high jet flow rate liquid 
surface breakup occurs and Sonin et al’s 10 condensation 
rate correlation does not apply. Therefore, the present 
data suggests that Sonin et al.’s 10 steady state condensa- 
tion rate correlation may be used to predict the tran- 
sient condensation rate in a mixing process if the 
instantaneous values of the jet subcooling and the aver- 
age turbulence intensity at the interface are used. 

Conclusion 

Experiments were conducted to investigate the mix- 
ing time and transient vapor condensation rate at the 
liquid-vapor interface in a thermally stratified liquid 
hydrogen tank mixed by an axial jet. Two mixing time 
definitions, based on the thermal equilibrium (t m ) th and 
pressure equilibrium (t m ) p , have been used. Under the 
conditions and parameter ranges in the present study, 
the following conclusions can be drawn: 

1. The Richardson number (Ri) defined in the pres- 
ent study is an appropriate parameter to characterize 
the effect of buoyancy in the determination of mixing 
times (t m ) th and (t m ) p . If the Richardson number is less 
than 0.5, the buoyancy effect is neglected. Therefore, the 
dimensionless mixing time for the high jet flow rates in 
a normal- gravity condition may be applied to the low 
jet flow rates in a zero-gravity environment if the effect 
of liquid-vapor interface configuration is neglected. 

2. For Ri < 0.5, the mixing times (t m ) th and (t m ) p 
can be predicted by Eqs. (12a) and (14a), respectively. 
These equations may also be applied to a space syatem 
in which the gravity effect is negligible. Buoyancy delays 
the mixing process. Equations (12b) and (14b) give the 
predictions of (t m ) th and (t m ) p , respectively, for Ri 
> 0.5 based on the system used the the present study. 

3. For a steady state condensation process with 

negligible buoyancy, the liquid subcooling (T s — T b ) can 


be expressed by Eq. (18) as a function of the jet subcool- 
ing (T g — Tj) and the average turbulence intensity (v g ) 
at the interlace. 

4. The steady state condensation rate correlation of 
Eq. (19) may be used to predict the transient condensa- 
tion rate in a mixing process if the instantaneous values 
of jet subcooling and average turbulence intensity at the 
interface are employed. However, if the liquid surface 
breakup occurs due to high jet flow rate, Eq. (19) may 
underpredict the vapor condensation rate. 
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TABLE 1.— THE INITIAL AND FINAL CONDITIONS OF THE TEST 


Test 

Liquid fill. 

Initial conditions 

Final 

Dome 

Jet flow 

series 

percent 




conditions, 

temperature, 

rate, 






P, 


9i’ 


by height 

by volume 

p, 


kPa 

K 

m 3 /hr 




kPa 

mmmrn 




A1 

76 

85 

235.1 

6.1 

43.4 

63.9 

0.568 

A2 

74 

83 

202.7 

6.1 

38.6 

67.2 

.772 

A3 

73 

81 

238.5 

6.2 

40.0 

75.6 

.681 

B1 

72 

80 

199.2 

5.5 

42.0 

77.8 

.409 

B2 

71 

78 

200.6 

5.4 

36.5 

82.2 

1.045 

B3 

70 

77 

187.5 

5.2 

36.5 

85.0 

1.635 

B4 

68 

75 

188.2 

5.2 

37.9 

88.9 

2.430 

Cl 

42 

44 

184.0 

4.9 

57.9 

129.4 

.409 

C2 

41 

42 

188.9 

5.0 

57.9 

137.2 

.681 

D1 

66 

73 

218.5 

6.2 

46.9 

91.1 

.454 

D2 

64 

71 

208.2 

5.9 

41.4 

96.1 

1.158 

D3 

63 

70 

197.8 

5.7 

42.0 

97.2 

1.863 


TABLE 2. — THE MIXING CHARACTERISTICS OF THE TEST 


Test 

series 

Qj, 

m 3 /hr 

Rej 

Ri 

H/D 

OiaJih 1 

min 

(Op* 

min 

IBM 

(O p 

A1 

0.568 

38,600 

1.68 

0.968 

3.73 

3.54 


6.18 

A2 

.772 


.77 

.943 

2.06 

1.60 


3.80 

A3 

.681 

47,600 

1.09 

.928 

2.37 

2.27 

5.15 

4.75 

B1 

.409 

27,900 

2.62 

.916 

6.60 

6.77 

8.67 

8.52 

B2 

1.045 

70,900 

.40 

.903 

1.33 

1.03 

m 

3.31 

B3 

1.635 

111,500 

.15 

.877 

.80 

.60 

1 

3.02 

B4 

2.430 

166,300 

.07 

.852 

.55 

.40 

4.45 

2.99 

Cl 

.409 

28,400 

.88 

.526 

2.80 

4.13 

4.86 

5.20 

C2 

.681 

46,900 

.31 

.513 

1.37 

1.80 

4.01 

3.77 


.454 

29,900 

2.09 

.826 

4.60 

4.47 

7.07 

6.24 

D2 

1.158 

77,200 

.30 

.801 


.97 

4.10 

3.45 

D3 

1.863 

124,500 

.11 

.789 


.56 

4.26 

3.21 
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Figure 1 .—The sketch for the test tank and the locations of tem- 
perature sensors and thermistors in the tank. 
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Tank Pressure, kPa Tank Pressure, kPa 




Figure 2.— The change of tank pressure as a function of time. 
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Time, Min. 

(d) Test series D. 
Figure 2. — Concluded. 
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Liquid Temperature, K Liquid Temperature, 


24 




Time, Min. 

(b) Test series C, Qj « 0.681 m 3 /hr. 

Figure 3.— The change of liquid temperature as a function of time. 
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Richardson Number, Ri 


Figure 4.— -The dimensionless mixing time (based on the 
thermal equilibrium) as a function of Richardson number. 



0 0.5 1 1.5 2 2.5 3 


Richardson Number, Ri 

Figure 5.— The dimensionless mixing time (based on 
the pressure equilibrium) as a function of Richardson 
number. 


IS 






0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Time, Min. Time, Min. 

(a) Tost series B. (b) Test series D. 

Figure 6. — The change of the Jet subcooling as a function of time. 
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Dimensionless Height, h/L Dimensionless Height, h/L 



(a) Test series B, Qj * 1 .045 m 3 /hr. 



Fluid Temperature, K 
(b) Test series B, Qj » 1 .635 m 3 /hr. 

Figure 7.—' The fluid temperature distribution as a function of height. 
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Vapor Mass, kg Vapor Mass, kg 



Time, Min. 

(a) Test series B, Qj = 1 .045 m 3 /hr. 



Time, Min. 


(b) Test series B, Qj = 1.635 m 3 /hr. 




Time, Min. 



Time, Min. 


(d) Test series D,Qj = 1.158 m 3 /hr. (e) Test series D, Qj = 1 .863 m 3 /hr. 

Figure 8.— The change of vapor mass as a function of time and its comparison with the prediction based on the condensation 
rate correlation equation. 
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